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streaming  Reynolds  number  of  approximately  240.  Within  the  laminar  regime,  the 
transition  from  “intermediate”  to  “large”  values  of  the  streaming  Reynolds  number 
occurs  at  approximately  130.  Heat  transfer  results  for  large  values  of  the  streaming 
Reynolds  number  in  the  laminar  regime  closely  match  the  present  theory  (less  than 
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I.  INTRODUCTION 


The  science  of  convective  heat  transfer  in  an  acoustic  field,  while  still  in  its  infancy, 
presents  many  new  and  exciting  possibilities  for  application  in  the  future.  Different 
experiments  and  theory  have  proven  that  under  the  correct  circumstances  there  are  certain 
desirable  heat  transfer  characteristics  affecting  an  object  which  is  immersed  in  a  strong 
acoustic  field.  A  complete  analysis  concerning  the  processes  and  effects  of  the  related  heat 
transfer  phenomena,  though,  is  lacking  and  desperately  needed. 

Although  thermoacoustics  have  already  been  applied  to  some  advanced  heat  transfer 
designs,  for  instance,  the  thermoacoustic  space  refrigerator  and  other  thermoacoustic 
cryocoolers  developed  at  the  Naval  Postgraduate  School,  there  has  yet  to  be  developed 
anything  that  can  compete  on  an  economic  level  with  what  is  currently  marketed  today.  The 
efficiencies  obtained  so  far  have  been  quite  low,  requiring  nearly  twice  the  power  of  a 
conventional  vapor  compression  refrigerator.  When  the  fundamentals  behind  the 
thermoacoustic  phenomenon  and  the  related  heat  transfer  characteristics  are  completely 
understood,  breakthroughs  can  occur  which  could  allow  industry  to  move  ahead  and  apply 
these  techniques  on  an  every  day  basis  toward  a  variety  of  common  uses. 

To  properly  model  and  control  the  parameters  which  impact  upon  the  heat  transfer 
behavior  in  a  thermoacoustic  engine,  it  would  be  advantageous  if  the  various  flow  regimes 
(e.g.,  turbulent  vs.  laminar)  in  the  engine  could  be  isolated  and  analyzed  in  detail.  Further 
information  in  this  regard  may  be  obtained  through  a  parametric  analysis  of  a  suitable  model 
problem  by  which  a  measure  of  the  importance  (or  rather  a  magnitude  of  the  effect  each 
parameter  has  on  the  process  as  a  whole)  can  be  determined.  This  is  an  important  element 
of  the  modeling  process  and  requires  study  since  by  understanding  the  impact  that  each 
individual  parameter  makes  upon  the  thermoacoustic  process  as  a  whole,  it  may  be  possible 
to  predict  the  changes  in  the  heat  transfer  characteristics  as  individual  components  are  varied. 

The  work  contained  within  provides  an  experimental  study  of  some  of  the  dominant 
heat  transfer  properties  of  a  particular  model  problem  that  may  be  encountered  in 
thermoacoustic  engines.  The  model  problem  chosen  is  one  of  convective  heat  transfer  from 
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a  cylinder  in  a  zero-mean  oscillatory  flow.  The  flow  is  representative  of  the  acoustic 
standing  wave  in  a  thermoacoustic  engine  whereas  the  cylinder  represents  a  tube  or  other 
component  that  may  be  present  in  such  an  engine. 

The  work  involves  a  correlation  of  experimental  heat  transfer  data  in  terms  of  a 
suitable  Nusselt  number  (Nu)  with  other  appropriate  dimensionless  parameters  in  the 
problem,  such  as  the  streaming  Reynolds  number  (Rs),  which  itself  is  a  function  of  length 
scales,  pressure  ratios  and  frequency  parameters  and,  of  course,  of  the  Prandtl  number  (Pr). 
Through  use  of  high  power  standing  resonant  acoustic  waves  in  a  cylindrical  chamber,  a 
high-intensity  internal  oscillatory  flow  is  established.  Under  these  conditions,  the  heat 
removal  rate  from  a  thin  cylindrical  heating  element  immersed  in  the  acoustic  signal  will 
supply  data  necessary  to  arrive  at  some  basic  conclusions  as  to  heat  transfer  phenomena 
occurring  around  a  cylinder.  This  work  focuses  upon  two  different  flow  regimes;  one  in 
which  laminar,  attached  flow  around  the  cylinder  at  large  values  of  the  streaming  Reynolds 
number  is  present,  and  the  second  in  which  vortex  shedding  and  other  instabilities  in  the  flow 
are  expected  to  occur  at  the  cylinder  surface.  The  resultant  experimental  data  additionally 
provides  guidelines  for  determining  when  the  flow  transitions  from  one  regime  to  the  other. 
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II.  BACKGROUND 


A.  HISTORICAL 

It  has  long  been  understood  that  a  large  temperature  gradient  along  the  length  of  a 
cylindrical  tube  can,  under  certain  suitable  circumstances,  spontaneously  excite  the  fluid  into 
oscillations  strong  enough  to  create  audible  sound.  Glass  blowers  provided  the  earliest 
accounts  of  this  acoustic  effect.  They  found  that  when  one  end  of  a  glass  tube  was  placed 
in  a  furnace,  the  temperature  difference  between  the  end  of  the  tube  in  the  furnace  and  the 
end  still  under  ambient  conditions  created  an  audible  tone  which  was  emitted  from  the  open 
end  of  the  tube.  Even  though  early  scientists  knew  of  this  effect,  it  was  merely  considered 
to  be  more  of  an  oddity  than  a  scientific  discovery  which  might  have  useful  implications. 
In  fact,  the  earliest  accounting  of  what  may  have  been  the  first  ever  thermoacoustic  engine, 
the  Sondhauss  tube  (Figure  1),  was  described  by  Sondhauss  (1850)  himself  as  the  “glowing 
glass  harmonica”.  Lord  Rayleigh  (1945)  gave  the  first  good  qualitative  analysis  of  the 
Sondhauss  tube  in  1 896  in  which  he  described  the  mechanism  behind  the  effect  and  posed 
the  theory  that  mechanical  work  could  be  obtained  from  the  vibrations,  or  oscillations,  which 
were  being  created  by  the  temperature  gradient  along  the  tube. 
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An  extension  of  research  into  Sondhauss  tubes  was  conducted  by  Carter  (1988)  and 
included  placing  a  stack  of  plates  at  an  appropriate  point  in  the  tube.  The  plates  included  hot 
heat  exchanger  strips  at  one  end  and  cold  heat  exchanger  strips  at  the  other.  These  improved 
the  effect  of  the  Sondhauss  tube,  and  inspired  another  scientist,  Feldman  (1966)  to  conduct 
similar  research  which  consequently  resulted  in  an  oscillator  which  produced  27  W  of 
acoustic  power  from  600  W  of  thermal  power.  All  of  this  occurred  in  the  1960s. 

It  wasn’t  until  much  later  that  it  was  postulated  that  the  effect  could  somehow  be 
reversed,  i.e.,  a  temperature  gradient  could  be  created  along  a  tube  using  a  powerful, 
resonating  acoustic  signal  as  the  driving  force.  It  is  only  more  recently  that  advances  in 
acoustic  technology  have  allowed  serious  research  toward  this  goal. 

B.  RECENT  WORK 

The  heat  transfer  effect  related  to  this  thermoacoustic  phenomenon  has  since  been 
used  in  the  creation  of  heat  pumps  and  has  been  explored  by  several  scientists.  The  results 
of  these  attempts  include  what  was  termed  a  pulse-tube  refrigerator  as  developed  by  Gifford 
and  Longsworth  (1966).  This  refrigerator  utilized  low  frequency,  high  amplitude  oscillations 
to  excite  the  gas  in  a  tube  and  create  a  cooling  effect  along  the  surface.  The  invention  of 
“modem”  thermoacoustic  refrigeration  occurred  in  the  early  1980s  at  Los  Alamos  National 
Laboratory.  It  was  in  essence  a  modification  to  the  work  that  Carter  did,  using  stacks  of 
plates  with  a  much  smaller  temperature  gradient.  Additional  engineering  developments  by 
others,  such  as  the  work  at  the  Naval  Postgraduate  School  with  a  thermoacoustic  refrigerator 
intended  for  use  on  the  space  shuttle,  led  to  increases  in  efficiency,  as  well  as  to  an  increase 
in  commercial  interest  and  development.  Currently,  there  are  major  projects  ongoing  in 
several  countries  as  outlined  by  Swift  (1995),  including  a  prototype  food  refrigerator  based 
upon  the  Naval  Postgraduate  School’s  work  being  built  in  the  Republic  of  South  Africa.  The 
Ford  Motor  Company  has  developed  its  own  version  of  a  thermoacoustic  refrigerator  while 
the  Tektronix  Corporation  is  working  towards  a  pulse-tube  type  of  refrigerator  to  be  used  for 
cooling  electronics  to  cryogenic  temperatures. 

All  of  these  attempts  are  especially  significant  given  the  current  stigma  surrounding 
the  environmentally  hazardous  use  of  CFC’s.  Even  though  thermoacoustic  refrigeration  is 
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still  not  as  efficient  as  that  of  current  energy  efficient  vapor  compression  models,  there  is  a 
growing  demand  for  something  to  take  their  place.  By  advancing  our  knowledge-base  in  this 
area,  and  further  incorporating  a  new  understanding  of  how  increases  to  the  efficiency  of 
thermoacoustic  designs  can  be  made,  that  it  may  well  be  possible  that  a  new  thermoacoustic 
revolution  is  in  our  future. 

C.  THERMOACOUSTIC  PROCESS 

The  basic  process  behind  a  thermoacoustic  engine  can  be  best  described  by  the  model 
in  Figure  2.  The  upper  portion  of  the  figure  shows  a  sound  chamber  with  an  acoustic  driver 
at  the  left  end  which  is  used  to  create  a  resonant,  standing  wave  in  the  chamber.  At  an 
appropriate  point  within  the  chamber,  a  thermoacoustic  stack  (Figure  3)  is  placed  with  a  heat 
exchanger  (Figure  4)  on  either  side  of  it.  The  flow  of  heat  is  from  right  to  left  in  the  figure. 
The  process  of  heat  transport  across  the  plate  is  illustrated  in  the  lower  portion  of  the  figure. 
The  fluid  within  the  chamber  will  oscillate  due  to  the  acoustic  wave,  traveling  from  a  point 
of  low  pressure  to  one  of  high  pressure,  gaining  and  losing  energy  during  each  half  cycle. 

For  instance,  a  parcel  of  gas  at  temperature  T0  at  low  pressure  moves  along  line  1, 
increasing  in  pressure  and  temperature  until  it  reaches  T++  (has  gained  two  units  of  heat). 
At  that  point,  it  loses  one  unit  of  heat  to  the  plate,  thereby  reducing  its  temperature  to  T+. 
As  the  parcel  of  gas  continues  through  the  second  half  of  the  oscillatory  cycle,  it  decreases 
in  pressure,  and  it  loses  two  unit  of  heat,  dropping  to  T-.  It  is  then  able  to  retrieve  a  single 
unit  of  heat  from  either  the  plate  at  the  right  hand  side  of  the  cycle,  or  from  a  heat  exchanger 
at  that  end.  In  effect,  a  “bucket-brigade”  of  little  parcels  of  gas  is  formed  as  the  heat 
transport  mechanism.  It  is  the  heat  exchangers  on  either  end  of  the  thermoacoustic  stack  and 
the  heat  transport  mechanisms  involved  with  them  that  this  experiment  intends  to  analyze. 

D.  HEAT  TRANSFER  IN  ACOUSTIC  FIELDS 

The  early  1950's  saw  an  increased  interest  towards  understanding  the  heat  transfer 
behavior  in  oscillatory  flows  and  an  earnest  effort  towards  understanding  the  possible 
benefits  thereof  began.  Richardson  (1967)  produced  the  first  significant  contribution  to  this 
field  by  providing  a  coherent  and  detailed  account  on  the  general  nature  of  heat  transfer  in 
oscillatory  flows.  He  gave  concise  documentation  on  how  sound  and  vibration  fields  had 
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made  an  impact  in  the  heat  transfer  field  up  to  that  point.  He  additionally  conducted  a 
preliminary  study  on  how  the  effect  can  be  related  to  that  of  a  more  traditional  fluid- 
mechanical  analysis.  Davidson  (1974)  further  expanded  upon  Richardson’s  work,  analyzing 
the  heat  transfer  behavior  of  cylinders  in  oscillatory  flow.  That  occurred  in  the  early  1970's, 
and  since  then,  very  little  work  on  the  subject  can  be  found  in  literature. 

Most  recently,  though,  Mozurkewich  (1995)  performed  heat  transfer  experiments  in 
an  acoustic  field  utilizing  a  transient  analysis  of  a  heated  wire.  His  results,  although 
informative  in  some  respects,  were  lacking  in  data  within  the  most  basic  flow  regimes  and 
his  conclusions  left  some  doubt  to  the  reader. 

When  first  attempting  to  analyze  the  process  behind  convective  heat  transfer  in  an 
acoustic  field,  it  is  often  simpler  to  think  of  the  heat  transfer  phenomenon  as  something  akin 
to  that  of  forced  convection  due  to  a  steady  mean  flow.  Initially,  that  there  is  a  separate 
power  source  placed  away  from  the  test  object  which  produces  a  disturbance  in  the  fluid 
medium  in  which  that  object  is  immersed.  In  forced  convection,  that  power  source  may  be 
considered  to  be  a  fan  or  a  pump  which  creates  a  pressure  gradient,  which  in  turn  causes  a 
steady  flow  of  fluid.  In  the  current  problem,  the  power  source  is  instead  an  acoustic  driver 
which  causes  an  oscillatory  (or  vibrational)  type  of  time-periodic  flow  around  the  object 
being  considered.  This  oscillatory  flow  has  a  zero-mean  and  results  in  no  net  through  flow. 
Before  analyzing  this  flow  for  heat  transfer  characteristics,  though,  it  is  necessary  to  first  note 
which  aspects  of  the  acoustic  field  dominate. 

A  resonant,  standing  wave  acoustic  field  excited  across  the  ends  of  a  closed 
cylindrical  chamber  has  very  distinct  properties.  Of  particular  interest  is  the  fact  that  when 
such  a  field  exists,  the  pressure  along  the  length  of  the  chamber  varies  sinusoidally  such  that 
a  point  of  maximum  pressure  occurs  at  the  rigid  end  termination  at  the  opposite  end  of  die 
chamber  from  which  the  acoustic  signal  is  being  generated.  In  addition,  depending  upon  the 
frequency  being  used,  there  may  be  more  than  one  zero-crossing,  or  pressure  minimum, 
along  the  length  of  the  chamber  (Figure  5). 
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At  resonance,  the  acoustic  velocity  is  out  of  phase  with  the  pressure.  For  instance, 
at  a  point  of  minimum  pressure,  a  pressure  node,  the  particle  oscillations  will  be  at  their  point 
of  maximum  velocity,  a  velocity  antinode.  The  reverse  then  also  holds  true,  (i.e.,  a  pressure 


Circular  Sound 


Velocity  Profile 
Pressure  Profile 


Figure  5.  Pressure  and  Velocity  Profiles  in  a  Circular  Tube  in  the 
Presence  of  an  Acoustic  Signal. 


antinode  can  also  be  designated  a  velocity  node). 

To  understand  the  reason  why  this  relationship  is  extremely  valuable  to  this  research, 
it  is  useful  to  again  refer  to  the  forced  convection  model.  In  general,  as  the  fluid  velocity  in 
a  forced  convection  application  increases  within  a  particular  regime,  so  does  the  heat  transfer 
coefficient.  The  case  is  essentially  the  same  for  heat  transfer  in  an  acoustic  field,  and  hence 
it  is  expected  that  the  heat  transfer  rate  would  also  increase  as  the  particle  velocity  increases. 
Thus,  for  experimental  purposes  it  is  of  prime  importance  that  the  test  object  be  placed  at  a 
velocity  anti-node  to  get  the  maximum  effectiveness  out  of  the  process  for  a  given  acoustic 
signal.  However,  in  contrast  to  the  forced  convective  mean  flow  case,  the  current  issue  with 
oscillatory  flow  is  considerably  more  complicated  due  to  the  wide  range  of  flow  parameters, 
and  hence  flow  patterns  and  heat  transport  mechanisms,  that  can  result. 
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II.  EXPERIMENT 


A.  INTENT 

As  with  any  other  advance  in  technology  or  new  scientific  discovery  with  which 
engineers  desire  to  predict  and  quantify  results  in  some  manner,  it  is  best  to  begin  by  first 
breaking  the  process  down  into  its  basic  component  parts.  By  completing  an  analysis  for  the 
simplest  version  of  the  model  in  question,  a  stepping  stone  will  be  established  upon  which 
the  analysis  of  more  complicated  scenarios  can  be  built.  It  stands  to  reason  that  this  is  how 
the  analysis  of  the  thermoacoustic  heat  transfer  process  should  also  begin. 

For  an  initial  starting  point,  the  analogy  once  again  to  forced  convection  heat  transfer 
is  used.  The  problem  of  an  isolated  cylinder  in  a  mean  cross  flow  is  well  documented  and 
understood,  and  has  in  turn  been  used  to  develop  correlations  for  flow  over  a  collection  of 
cylinders,  such  as  a  tube  bank,  a  more  practical  application  as  is  evident  from  any  basic  heat 
transfer  textbook.  This  is  the  motivation  for  a  study  of  the  behavior  of  simple  shapes  such 
as  a  cylinder  placed  in  an  acoustic  field.  It  is  hoped  that  with  this  knowledge  for  an  isolated 
cylinder,  it  would  be  possible  to  extend  the  solution  to  other  models. 

It  is  only  necessary  then  to  concentrate  on  the  evaluation  of  the  acoustic  signal  itself 
in  terms  of  its  many  different  parameters.  But  before  this  analysis  may  begin,  the  type  of 
acoustic  flow  around  the  cylinder  must  be  established.  To  allow  corroboration  with 
established  theory,  the  flow  pattern  initially  desired  is  that  of  basic  laminar,  incompressible 
flow  where  well  understood  streaming  patterns  are  the  principal  forms  of  flow  and  heat 
transport.  This  meets  the  requirement  for  maintaining  the  most  simplified  version  of  flow 
for  the  analysis. 

The  geometry  of  how  the  test  cylinder  is  placed  with  reference  to  the  acoustic  signal 
is  also  of  utmost  importance.  A  theoretically  perfect  scenario  would  have  the  test  cylinder 
situated  normal  to  a  unidirectional  sound  field  with  no  interference  from  the  surroundings. 
Such  a  situation  cannot  be  exactly  duplicated,  though,  due  to  the  nature  of  acoustic  waves 
to  spread  and  travel  in  all  directions.  In  order  to  restrict  the  acoustic  signal  to  only  one 
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direction,  it  is  proposed  that  the  sound  field  be  set  up  in  a  cylindrical,  resonant  sound 
chamber  so  that  the  analysis  is  properly  limited  to  axial  wave  modes  only. 

B.  FUNDAMENTAL  IDEAS 

Unlike  the  well  known  dependence  on  the  Reynolds  and  Prandtl  numbers  found  in 
conventional  cross  flow  over  bluff  bodies,  the  issue  of  heat  transfer  in  the  presence  of  an 
acoustic  field  is  significantly  more  complicated  due  to  the  presence  of  a  multitude  of 
competing  length  (or  time,  or  velocity)  scales.  The  ways  in  which  these  length  scales  may 
be  ordered  are  many  and  lead  to  numerous  distinct  parameter  regimes  with  quite  drastically 
different  flow  properties,  and  hence  heat  transfer  properties.  In  order  to  closely  examine  the 
properties  of  heat  transfer  in  oscillatory  flows  then,  it  is  first  necessary  to  enumerate  some 
of  the  different  parameters  and  variables  involved.  By  establishing  and  maintaining  a  set  of 
criteria  surrounding  these  parameters,  and  by  modeling  the  test  apparatus  to  conform  to  them, 
the  job  of  evaluating  different  regimes  of  flow  and  drawing  correlations  from  the  data 
obtained  will  become  possible. 

This  section  lists  the  most  important  of  these  parameters  and  provides  reasoning  for 
the  choices  involved,  keeping  in  mind  the  desire  to  maintain  the  most  basic  and  core  set  of 
conditions  that  flow  around  the  test  cylinder  be  laminar,  incompressible  and  attached,  and 
that  only  the  steady-state  solution  be  considered. 

1.  Criterion  A 

Following  the  first  assumption  that  the  acoustic  streaming  flow  around  the  test 
cylinder  be  incompressible,  two  different  length  scales  must  be  specified.  The  first  of  these 
requires  that  the  relative  size  of  the  test  cylinder  be  very  small  compared  to  the  radian 
wavelength  of  the  acoustic  field,  where  the  radian  wavelength  is  defined  as 


and  the  characteristic  length  scale  of  the  test  cylinder  is  chosen  as  the  radius,  a.  The  radian 
wavelength  can  be  related  to  the  frequency  by: 
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(2) 


X  -  —  -  s!L  -  c 

2 tz  2%  o) 

Now,  by  asserting  that a«A,  and  designating  the  ratio  between  the  two  as  %,  it  follows  that: 


(3) 


and  finally  that 


(4) 


The  criterion  %  «  1  ensures  that  radiation  effects  due  to  the  acoustic  streaming  are 
negligible,  as  presented  by  Lighthill  (1963)  and  indicates  that  it  is  only  the  local  acoustic 
field  conditions  that  are  of  importance. 

2.  Criterion  B 

The  second  criterion  which  is  required  to  support  the  assumption  that  flow  be 

incompressible  is  derived  from  the  relationship  between  displacement  amplitude  of  particle 

—  "T 

oscillation  in  the  sound  field,  A ,  and  the  cylinder  radius,  a.  The  ratio  between  the  two,  — , 

a 

will  be  designated  as  the  amplitude  parameter,  e,  and  dictates  whether  or  not  separation  will 
occur. 

When  the  amplitude  parameter  is  very  small, 

e  =  —  «  1  (5) 

a 


the  particles  in  the  sound  field  move  a  very  short  distance  along  the  cylinder  wall  before 
reversing  their  direction.  This  ensures  that  the  flow  remains  attached,  with  little  chance  of 
separation  occurring,  and  hence  the  flow  will  remain  laminar  at  all  times.  It  can  also  be 
noted  that  eis  directly  proportional  to  the  pressure  ratio  P0  /  Pm  and  can  therefore  take  on 
much  larger  values  for  strong  acoustic  fields.  This  is  accomplished  by  observing  that  the 
displacement  amplitude  of  a  particle  in  the  flow  is  directly  related  to  its  velocity  and  that  the 
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particle  velocity  is  in  turn  related  to  the  pressure  ratio  by  the  following  (in  a  plane  standing 
sound  field) 


(O 


and 


cPq 

Y K 


(6) 


(7) 


therefore 


— «  1 
ypm 


(8) 


Yet  another  form  of  this  parameter  often  used  in  the  literature  is  called  the  Keulegan 
Carpenter  number  and  is  defined  as  KC-UJ2af  and  can  additionally  be  expressed  as 


KC  -  7U€  (9) 

Of  importance  is  the  fact  that  the  product  of  the  parameters  defined  in  the  first  two  criteria 
(A  and  B)  is  the  flow  Mach  number,  which  can  be  defined  as 

U0 

M=xe=—  (10) 

c 

When  criteria  A  and  B  are  satisfied,  M«l,  and  this  in  turn  is  the  second  condition  which 
satisfies  the  assumption  of  incompressible  flow. 

3.  Criterion  C 

The  Stokes  boundary  layer  thickness  S  is  related  to  the  kinematic  viscosity  and  the 
radian  frequency  of  oscillations  by 
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and  is  the  well  known  length  scale  which  is  a  measure  of  the  extent  of  viscous  effects  in  an 
oscillatory  flow.  A  frequency  parameter  A2  can  be  defined  as  follows 

A2  =  (j)2  =  —  (12) 

6  v 

For  the  case  when  A2»l,  the  Stokes  shear  layer  is  confined  to  a  narrow  region  and  the 
acoustic  streaming  effect  appears  as  a  slip  velocity  along  the  cylinder  surface.  Utilizing  the 
knowledge  that  the  boundary  layer  thickness  is  on  the  order  of  105  and  imposing  the 
condition  (somewhat  arbitrarily)  that 


a 

106 


>  4 


(13) 


it  follows  that 


A2  >  1600  (14) 

is  a  good  criterion  to  ensure  “large”  values  of  the  frequency  parameter. 

The  frequency  parameter  may  also  be  often  found  in  the  literature  in  the  form  of 
P  =(2 a)2f/\ ,  and  can  be  expressed  as 


P  =  (-)A2 

7t 


(15) 


4.  Criterion  D 

When  criteria  A  -  C  are  satisfied,  the  acoustic  streaming  velocity  is  of  magnitude 
O(eU0).  A  streaming  Reynolds  number,  Rs,  can  then  be  defined  as 


(eU0)a 


(16) 
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Through  substitution  for  e  and  U0 


which  yields 


R. 


U£ 

a(j)  v 


cov 


— (— )2 
WV  Y Pm 


Rs  =  e2A2 


(17) 


(18) 


This  streaming  Reynolds  number  becomes  the  driving  factor  in  determining  what  will  be  the 
primary  mode  of  heat  transport  within  the  region  due  to  the  acoustic  streaming  flow.  Stuart 
(1966)  demonstrated  that  when  Rs«l,  a  Stokes  flow  becomes  prevalent  in  the  outer  region 
while  a  boundary  layer  flow  is  predominant  when  this  parameter  takes  on  values  much 
greater  than  one.  In  order  to  ensure  that  forced  convective  heat  transfer  is  dominant  then,  we 
impose  the  following  constraint 

Rs  »  1  (19) 


5.  Criterion  E 

It  was  first  observed  by  Honji  (1981)  that  flow  around  a  cylinder  will  become 
centrifugally  unstable  and  separate  into  vortices  as  the  amplitude  of  particle  oscillation 
increases.  This  instability  occurs  in  the  Stokes  layer  where  the  flow  is  parallel  to  the  direction 
of  particle  oscillation.  This  was  confirmed  by  Hall  (1984)  who  conducted  a  linear  stability 
analysis  on  the  unsteady  boundary  layer  in  the  high-frequency  limit.  He  found  that  a  critical 
value  of  the  Reynolds  streaming  number  exists  for  which  instabilities  begin  to  form,  namely 
when  Rs  becomes  greater  than  4.24 A.  Recent  work  in  the  area  by  Sarpkaya  (1986)  provides 
further  explanation  for  this  phenomenon.  Since  vortex  shedding  can  make  a  large  impact  on 
the  heat  transport  from  the  cylinder,  we  will  maintain  the  following  criterion,  which  can  be 
expressed  in  two  different  ways  as 


Rs  <  4.24A 


(20) 
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or,  alternately  as 


v/A 

An  example  of  what  these  vortices  may  look  like  is  shown  in  Figure  6. 

6.  Criterion  F 

In  order  to  maintain  the  condition  in  which  there  is  minimal  influence  from  the 
buoyancy  effects  of  natural  convection  as  compared  to  the  forced  convective  heat  transfer 
due  to  the  acoustic  streaming  effects,  a  dimensional  analysis  of  the  governing  equations 
produces  the  following  requirement  for  the  Grashof  number 


(22) 


For  the  case  when  Gr  /  R*  ~  1  or  greater,  buoyancy  effects  must  be  taken  into  account  and 
any  heat  transfer  correlations  developed  will  have  to  be  modified  accordingly. 


Figure  6.  Vortices  Due  to  Oscillatory  Flow  Around  a 
Cylinder.  (Sarpkaya,  1986) 


7.  Criterion  G 

This  final  criterion  stems  from  the  desire  to  maintain  only  a  plane  standing  sound 
field  within  the  test  chamber.  This  means  that  only  the  axial  wave  components  excited  by 
the  resonant  acoustic  signal  be  present,  without  interference  from  transverse  modes  such  as 
the  azimuthal  modes  and  radial  modes.  It  may  be  shown  from  theory  that  there  is  a  certain 
cut-off  frequency  developed  from  a  solution  of  the  wave  equation,  below  which  the 
transverse  modes  will  not  be  present.  For  the  cylindrical  waveguide  geometry  being  used, 
the  solution  for  the  above  condition  is  obtained  in  the  form  of  appropriate  roots  of  Bessel 
functions.  If  L  is  the  length  of  the  test  chamber,  the  longitudinal  frequency  can  be  expressed 
as 


(23) 


where  l  is  the  1st,  2nd, ... ,  mode  number.  The  transverse  frequency  mode  is  expressed  by 


£  1 

f  =  a 

J  mn  mi 


tzD 


(24) 


where  D  is  the  test  chamber  diameter  and  anm  represents  the  eigenvalues  obtained  from  roots 
of  the  Bessel  functions.  Transverse  modes  will  be  present  if  /  =  ,  so  it  is  desired  to 

maintain/;  well  below /^  and/;  now  represents  the  maximum  frequency  possible  while  still 
maintaining  this  criterion.  By  substituting  for  both  frequencies,  the  condition  becomes 


Jc 

2 L 


< 


a„ 


c 

tzD 


which  is  reduced  to 


/ 


2(A) 


O' 


■K 


By  introducing  a  new  parameter  called  the  aspect  ratio, 


(25) 


(26) 
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(27) 


and  finding  the  smallest  root  of  this  Bessel  function, 

/ 

— min  =  0.586 
n 


(28) 


the  criterion  now  becomes 


—  <  0.586 
2 Z 


(29) 


which  can  be  rearranged  to  finally  get 


where 


Z  = 


D 


>  0.85  L 


(30) 


L 


f 

J  m 


(31) 


This  gives  a  relationship  between  the  geometry  and  the  maximum  frequency.  But  since  the 
maximum  frequency  is  already  defined  using  criteria  A  -  F,  criterion  G  in  fact  gives  us  the 
maximum  diameter  that  can  be  used,  or  if  the  diameter  is  also  given,  it  determines  the  length 
of  the  chamber  instead. 

8.  Basic  Flow  Description 

When  a  cylinder  is  immersed  in  a  standing  acoustic  field  and  all  of  the  previous 
criteria  have  been  met,  particle  oscillations  will  initiate  the  most  basic  form  of  an  acoustic 
streaming  flow.  This  steady  flow  (Figure  7)  is  symmetrical  about  two  axis,  is  circular  in 
nature,  and  includes  a  well  defined  boundary  layer.  The  boundary  layer  is  actually  quite  small 
and  is  greatly  exaggerated  in  the  figure  for  clarity.  Although  this  is  not  the  only  flow  which 
is  present,  it  does  represent  the  flow  pattern  that  has  the  biggest  impact  as  a  heat  transport 
mechanism  when  these  criteria  are  satisfied. 
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Figure  7.  Outer  Acoustic  Streaming  Flow  (also  known  as  Boundary  Layer  Flow) 


As  discussed  earlier,  this  fluid  flow  regime  is  the  main  focus  of  this  investigation  in 
which  at  first  we  consider  only  laminar,  unseparated  flow  around  the  test  cylinder.  One  of 
the  goals  of  this  experimental  study  is  to  determine  the  range  of  values  for  the  previously 
discussed  parameters  where  one  flow  regime  transitions  to  another  and  relate  them  in  terms 
of  the  streaming  Reynolds  number.  It  is  believed  that  such  transitions  will  be  reflected  in  the 
heat  transfer  behavior. 

In  order  to  define  a  specific  range  for  laminar  flow,  though,  these  criteria  needed  to 
refinement  by  asking  the  questions;  how  small  is  “much  less  than  one”  or  how  big  is  “much 
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greater  than  one”?  It  was  determined  that  Criterion  A  in  Eq.  4  would  be  met  by  picking 
X  <  0.1  and  Criterion  B  in  Eq.  5  would  be  met  by  choosing  e  <  0.3.  This  ensures 
incompressible  flow  around  the  test  cylinder  and  that  flow  remains  attached.  Criterion  C 
becomes  A2  >  1600,  confining  the  Stokes  shear  layer  to  a  narrow  region.  The  above  ranges 
have  been  fixed  based  on  prior  experience  and  preliminary  experimentation,  but  are  not  by 
any  means  intended  to  be  “hard  and  fast”.  They  may  indeed  be  modified  if  necessary. 
Criteria  D  and  E  require  verification  after  each  experimental  run. 

C.  APPARATUS 

As  previously  discussed,  the  experiment  consists  of  a  heated  test  cylinder  placed 
normal  to  a  simple  acoustic  standing  wave  within  a  resonating  sound  chamber  (Figure  8). 
This  general  description  can  be  broken  down  into  three  major  components;  the  test  cylinder, 
the  sound  chamber  and  the  acoustic  electronics  package,  which  are  further  described  below. 


Figure  8.  Experimental  Test  Apparatus. 
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1.  Test  Cylinder 

It  was  originally  proposed  that  a  Watlow  stainless  steel  cartridge  heater  be  used  as 
the  test  cylinder  for  the  experiments.  They  are  available  in  varying  diameters  and 
additionally  feature  an  imbedded  type  “J”  thermocouple  placed  at  the  midpoint  of  its  length. 
This  presented  problems,  however,  since  the  heating  along  the  length  of  the  cartridge  heater 
was  uneven  and  the  relatively  large  thermal  resistance  of  the  stainless  steel  produced  large 
variances  in  surface  temperature  along  the  length.  Since  constant  surface  temperature  is  a 
feature  which  is  very  important  to  the  analysis  of  heat  transfer  characteristics,  this  was 
unacceptable.  Therefore,  a  copper  sheath  was  designed  to  fit  over  the  cartridge  heater.  It 
was  assumed  that  the  large  thermal  conductivity  of  copper  would  even  out  the  axial  surface 
temperature  gradient.  In  addition,  silicon  oil  was  liberally  applied  to  the  inside  of  the  copper 
sheath  prior  to  insertion  of  the  cartridge  heater  before  each  experimental  run  (Figure  9).  This 
provides  better  thermal  contact  in  the  narrow  annular  gap  (approximately  2  mm)  between  the 
sheath  and  heater,  preventing  air  gaps  which  can  cause  local  temperature  discontinuities  at 
the  surface.  This  arrangement  (as  was  later  verified.  Appendix  B)  provides  for  an  axial 
temperature  variation  of  less  than  0.5°C  from  one  end  of  the  heater  to  the  other  for  the 
temperature  range  of  interest  in  this  experiment.  A  picture  of  the  test  cylinder  is  shown  in 
Figure  10. 
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As  was  stated  earlier,  the  surface  temperature  of  the  test  cylinder  is  a  required  datum 
point  for  the  analysis.  Since  the  only  provision  for  temperature  measurement  is  at  the 
thermocouple  placed  at  the  center  of  the  test  cylinder,  the  equivalent  resistance  of  the  steel/ 
oil/copper  circuit  is  required  in  order  to  deduce  the  surface  temperature  from  the  cartridge 
heater  center  temperature  as  measured  by  the  imbedded  thermocouple.  Appendix  B  gives 
a  detailed  analysis  of  the  derivation  of  the  resistance,  which  is  approximately  1.019  K/W. 

The  cartridge  heater  receives  its  power  from  a  Kikusui  Model  PAR  160 A  regulated 
DC  power  supply.  It  provides  power  control  measurement  down  to  0.01  amps  and  0.01 
volts.  Thermocouple  measurements  are  provided  by  a  Keithley  Model  740  scanning 
thermometer  system.  Calibration  data  for  all  thermocouples  and  the  thermocouple  reader  is 
provided  in  Appendix  B. 

2.  Sound  Chamber 

The  purpose  of  the  sound  chamber  is  to  provide  an  environment  through  which 
acoustic  signals  of  various  frequencies  can  be  used  to  excite  a  resonant  standing  wave.  In 
order  to  accomplish  this  for  all  frequencies  which  may  be  used  during  the  experiment,  a 
resonant  chamber  that  would  be  adjustable  in  length  was  highly  desired.  Figure  1 1  shows 
the  final  configuration  of  the  test  chamber  while  Figure  12  shows  a  photo  of  the  test 
apparatus. 


Figure  11.  Sound  Chamber  Assembly. 
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Figure  12.  Photo  of  Test  Apparatus.  Acoustic  Driver  is  on  to  the  right. 


The  test  chamber  itself  is  a  plexiglass  tube  approximately  two  meters  long  which  is 
mounted  firmly  to  a  plexiglass  base  plate.  The  acoustic  driver  is  mounted  at  one  end  and 
provides  the  source  of  the  acoustic  signal.  The  opposite  end  has  a  movable  flat  end  plate 
which  has  an  o-ring  seal  to  help  isolate  the  acoustic  wave  in  the  chamber.  The  end  plate 
serves  as  the  rigid  end  termination  while  the  o-ring  provides  good  sound  confinement,  as 
well  as  good  stability  for  the  end  face  so  that  it  doesn’t  become  offset  to  either  side  as  its 
position  is  varied  along  the  length  of  the  enclosing  resonant  chamber. 

A  hole  through  this  end  plate  provides  access  for  a  pressure  transducer  which  is  used 
to  help  deduce  the  pressure  ratio,  and  hence  the  sound  pressure  level  for  that  configuration 
of  frequency  and  input  power.  The  pressure  transducer  used  is  an  Endevco  model  851  OB-5 
which  has  an  output  in  millivolts  and  has  a  pressure  sensitivity  of  50.89  mV/psi.  This  is 
connected  to  a  preamplifier  with  a  gain  of  100  and  provides  output  to  both  an  oscilloscope 
and  a  Hewlett-Packard  model  34401 A  multimeter.  Measurement  of  the  voltage  output  from 
the  multimeter  is  important  to  determine  when  resonance  has  occurred  within  the  chamber 
for  as  the  frequency  is  varied,  the  output  voltage  from  the  microphone  decreases  on  either 
side  of  the  resonant  operating  point. 

The  oscilloscope  provides  a  visual  representation  of  the  time  trace  of  the  acoustic 
signal  at  the  end  of  the  sound  chamber,  corresponding  to  a  pressure  antinode,  and  is  used  to 
ensure  that  the  signal  remains  sinusoidal  throughout  the  experimental  range  of  powers  and 
frequencies.  During  the  initial  stages  of  the  experiment,  the  oscilloscope  allowed  for  the 
discovery  of  interference  patterns  in  the  sinusoidal  waveform  as  caused  by  higher  order 
harmonics  at  high  SPLs  (>  155  dB).  In  order  to  limit  the  interference  that  was  present,  the 
use  of  “sleeves”  within  the  chamber  was  recommended  to  detune  the  resonant  mechanism. 
This  would  prevent  the  harmonics  from  being  integral  multiples  of  each  other  and  thereby 
prevent  them  from  reinforcing  each  other  to  form  interference  patterns.  By  placing  a  sleeve 
at  an  appropriate  spot  in  the  chamber,  some  of  the  high  frequency  harmonics  leading  to 
interference  could  be  eliminated,  allowing  for  even  higher  SPLs  to  be  achieved  before 
interference  occurred. 
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3.  Acoustic  Electronics  Package 

Of  great  importance  to  the  experiment  is  the  ability  to  generate  a  nearly  pure 
sinusoidal  waveform  at  varying  frequencies  and  high  power  ranges.  As  the  strength  of  the 
signal  generated  increases,  the  effect  of  the  flow  around  the  test  cylinder  becomes  more 
pronounced,  enhancing  the  heat  transfer  characteristics  of  the  system. 

The  acoustic  signal  being  generated  at  the  driver  end  of  the  sound  chamber  is 
provided  by  a  Hewlett-Packard  model  33120A  arbitrary  waveform  generator.  It  sends  a 
sinusoidal  waveform  at  the  proper  frequency  through  a  Techron  model  7540  power  supply 
amplifier  to  a  JBL  model  2490H  acoustic  compression  driver. 

A  more  detailed  review  of  each  item  of  equipment  used  in  the  system  is  provided  in 
Appendix  D. 

D.  EXPERIMENTAL  METHOD 

Prior  to  gathering  data  for  this  experiment,  it  was  first  necessary  to  develop  a 
coherent  plan  with  which  to  approach  the  problem.  The  first  step  required  was  to  find 
specific  frequencies  at  which  resonance  occurred  within  the  chamber,  and  which  would  also 
provide  a  velocity  anti-node  at  the  position  of  the  test  cylinder.  By  looking  at  the  geometry 
of  the  problem  (Figure  13),  it  became  obvious  that  these  limiting  factors  could  be  met  by  a 
combination  of  adjusting  the  end  plate  distance  from  the  test  cylinder,  as  well  as  the  signal 
frequency,  so  that  the  length  L  from  the  heated  cylinder  to  the  end  plate  termination  was  an 
odd  multiple  of  A/4.  This  could  be  further  refined  to  state  that  the  value  of  4Lf/c  needs  to  be 
an  odd  integer.  When  this  condition  was  satisfied,  the  requirement  of  a  velocity  anti-node 
occurring  at  the  cylinder  location  was  met. 

The  next  thing  needed  was  an  estimate  of  the  maximum  pressure  ratio,  and  therefore 
the  sound  pressure  level  that  could  be  obtained.  This  was  achieved  by  increasing  the 
amplitude  of  the  input  signal  waveform  until  the  output  waveform  on  the  oscilloscope  began 
showing  traces  of  interference  or  other  disturbances.  When  the  maximum  input  amplitude 
was  obtained,  the  multimeter  output  was  recorded. 
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where  n  =  5 
for  this  case 

x  — 


Figure  13.  Proper  Geometry  so  that  Maximum  Velocity  occurs  at  the  Test 
Cylinder. 


In  order  to  convert  the  output  voltage  to  an  actual  pressure  ratio  and  sound  pressure 
level,  it  was  necessary  to  first  understand  what  form  the  multimeter  uses  to  present  the 
output.  The  multimeter  gives  the  voltage  output  in  terms  of  the  true  RMS  value  of  the 
sinusoidal  signal,  or  rather  V^yjl,  where  Vo  represents  the  voltage  amplitude  from  zero  to 
peak  of  the  sinusoidal  signal.  Recall  that  the  true  RMS  value  for  the  pressure  can  also  be 
expressed  as  PJ\f2.  Then,  the  following  relationship  holds  true,  that 

Pjfi  =  (32) 

where  S  is  the  sensitivity  of  the  pressure  transducer  in  mV/Pa. 

The  sensitivity  of  the  pressure  transducer  after  passing  through  the  100  gain  setting 
of  the  preamplifier  was  converted  and  expressed  as  0.74  mV/Pa.  By  substituting  into  the 
above  equation,  the  pressure  amplitude  becomes 
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P  -  vJL 

0  0.523 mV/Pa 

where  it  is  noted  again  that  V0  I  \fl  is  the  multimeter  output. 
From  this,  a  pressure  ratio  can  be  defined  as 


(33) 


PR  = 


(34) 


where  Pm  is  the  mean  ambient  pressure  of  101  kPa.  The  pressure  ratio  is  typically  expressed 
in  terms  of  a  percentage. 

The  sound  pressure  level  is  defined  as  the  logarithmically  scaled  ratio  of  the  RMS 
pressure  and  a  predetermined  reference  pressure,  Pref,  where  Pref  is  chosen  to  be  20  pPa  for 
gases  (by  convention).  This  gives 


SPL  =  20 


log 


10 


ref 


(35) 


In  order  to  obtain  an  idea  of  how  strong  the  acoustic  signal  is  during  the  experiments.  Table 
1  gives  the  sound  pressure  level  for  various  activities 


Activity 

Sound  Pressure  Level 

Pressure  Ratio 

(dB) 

(%) 

Normal  Conversation 

60 

<0.001 

Jet  Airplane  at  Take-off 

90 

<0.01 

Pain  Threshold 

120 

0.28 

Minimum  Experimental  Level 

150.6 

0.94 

Maximum  Experimental  Level 

161.2 

3.20 

Table  1.  Samples  of  SPL  and  PR  for  Comparison. 


After  finding  the  pressure  ratio,  and  therefore  the  SPL,  the  actual  data  extraction 
phase  of  the  experiment  could  be  initiated.  It  was  preceded  by  with  an  understanding  of  what 
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type  of  data  was  needed  for  analysis.  A  crucial  element  of  this  experiment  is  the  need  to 
obtain  a  broad  base  of  data  with  which  to  incorporate  the  results.  In  order  to  obtain  this, 
experiments  were  run  at  several  different  pressure  ratios  for  each  resonant  frequency  found, 
starting  from  as  low  as  0.9%  and  building  up  to  the  maximum  pressure  ratio  by  increments 
of  0.1%.  In  order  to  obtain  a  spread  of  data  points  at  each  pressure  ratio  evaluated,  the  test 
cylinder  was  heated  to  approximately  8,  12  and  16  degrees  above  the  ambient  temperature. 
This  kept  the  power  requirements  low  and  reduced  the  amount  of  thermal  input  into  the 
ambient  air  within  the  sound  chamber.  This  latter  point  was  necessary  as  the  ambient 
temperature  could  rise  as  much  as  0.3  degrees  during  a  single  experimental  trial.  In  order  to 
ensure  reproducibility  and  reduce  the  effect  of  anomalistic  behavior,  three  runs  at  each 
temperature  point  were  conducted  to  ensure  consistency  of  data.  The  selected  temperatures 
were  only  guidelines  and  were  not  meant  to  be  hard  set  points  for  the  experiment.  Instead, 
they  were  treated  as  aim  points  with  an  acceptable  range  of  ±  1  degree.  Therefore,  in  order 
to  produce  an  even  broader  spread  of  data,  the  power  input  to  the  test  cylinder  was  varied 
slightly  for  each  of  the  three  trials  at  each  specific  temperature  point. 

Once  the  selected  frequencies  and  pressure  ratios  had  been  determined  and  a  suitable 
starting  pressure  ratio  and  temperature  had  been  obtained,  the  experimental  process  was 
initiated.  In  order  to  obtain  the  selected  pressure  ratio  for  a  particular  set  of  runs,  it  was 
necessary  to  get  an  idea  of  the  settings  required  for  each  specific  piece  of  equipment.  This 
was  done  by  selecting  the  appropriate  frequency  on  the  waveform  generator  and  modifying 
the  power  amplification  on  both  the  waveform  generator  and  the  power  amplifier  until  the 
appropriate  pressure  ratio  was  obtained  from  the  multimeter.  The  frequency  was  then 
adjusted  to  fine  tune  the  resonance.  A  check  of  the  oscilloscope  at  this  point  ensured  that  the 
signal  being  generated  was  of  the  right  waveform  and  that  interference  was  not  occurring. 
Then  the  power  amplifier  to  the  acoustic  driver  was  turned  down  to  zero  after  noting  the 
level  at  which  it  was  set.  This  allowed  for  obtaining  the  correct  pressure  ratio  in  a  quick 
manner  by  simply  turning  the  power  amplifier  up  to  the  previously  noted  value. 

It  was  necessary  at  this  point  to  ensure  that  the  test  cylinder  was  properly  prepared. 
This  entailed  introducing  approximately  ten  drops  of  silicon  oil  into  the  copper  sheath  and 
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inserting  the  cartridge  heater.  The  cartridge  heater  would  be  completely  immersed  in  the 
silicon  oil  after  being  fully  inserted.  After  a  period  of  time,  there  would  be  some  loss  of 
silicon  oil  due  to  the  wick  action  of  the  thermocouple  and  power  leads  emerging  from  the 
top.  This  was  insignificant  during  the  runs  required  for  a  single  pressure  ratio  and  caused  a 
negligible  change  in  the  center-to-surface  resistance  as  noted  in  Appendix  B,  but  it  became 
good  practice  to  add  some  of  drops  of  oil  each  time  a  new  set  of  runs  at  a  different  pressure 
ratio  were  to  be  taken. 

The  test  cylinder  was  then  inserted  into  the  sound  chamber,  making  sure  that  the 
bottom  of  the  cylinder  was  resting  in  a  shallow  indentation  specifically  machined  into  the 
inside  face  of  the  chamber.  Power  to  the  cartridge  heater  was  then  turned  on  and  set  so  that 
a  steady  state  temperature  of  approximately  eight  degrees  above  ambient  would  occur  when 
the  acoustic  signal  was  present.  This  became  more  of  an  art  form  and  required  familiarity 
with  the  system  to  accomplish  it  with  any  degree  of  accuracy. 

Since  the  power  to  the  acoustic  driver  at  this  point  was  still  at  zero,  the  power  to  the 
cartridge  heater  would  drive  the  temperature  of  the  test  cylinder  past  the  projected  steady 
state  temperature.  As  it  approached  the  projected  temperature,  though,  the  power  to  the 
acoustic  driver  was  increased  to  the  previously  noted  set  point.  This  provided  the  fluid  flow 
at  the  predetermined  pressure  ratio,  in  effect  beginning  to  cool  the  test  cylinder  until  it 
reached  a  steady  state  condition.  At  this  point,  the  rate  of  energy  input  to  the  cylinder  was 
equal  to  the  rate  of  energy  being  convected  away  from  the  cylinder.  By  monitoring  the 
interior  thermocouple  temperature,  it  was  easy  to  see  when  the  lowest  temperature  was 
reached.  When  the  temperature  began  to  rise  once  again,  it  was  determined  that  the  steady 
state  condition  had  been  reached  and  that  the  resultant  temperature  rise  being  witnessed  was 
due  only  to  the  test  cylinder  transferring  heat  into  the  surrounding  fluid  medium,  thus  raising 
the  overall  ambient  temperature. 

Two  situations  other  than  the  ideal  one  presented  above  occurred  frequently  due  to 
the  coarse  means  of  trying  to  arrive  at  the  desired  cylinder  temperature.  If  upon  engaging  the 
power  to  the  acoustic  driver  the  temperature  of  the  test  cylinder  did  not  decrease  but 
continued  to  increase  at  a  slower  rate  instead,  then  the  power  being  supplied  to  the  cartridge 
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heater  was  deemed  too  high  and  the  voltage  reduced  until  a  decrease  in  temperature  was 
witnessed.  If  upon  engaging  power  to  the  acoustic  driver  the  temperature  of  the  test  cylinder 
were  to  continue  decreasing  lower  than  the  desired  temperature  range,  it  was  an  easy  matter 
to  increase  the  voltage  to  the  cartridge  heater  to  increase  the  steady  state  temperature 
solution. 

When  the  steady  state  solution  point  was  reached,  the  frequency  and  microphone 
voltage  were  noted,  as  was  the  voltage  and  current  supplied  to  the  cartridge  heater  and  the 
temperature  of  the  thermocouple  in  the  cartridge  heater.  The  power  to  the  heater  and  the 
power  to  the  acoustic  driver  were  then  simultaneously  turned  off  and  the  test  cylinder  was 
quickly  removed  from  the  sound  chamber.  A  second  thermocouple  was  then  introduced 
through  the  access  hole  in  the  sound  chamber  (from  which  the  test  cylinder  had  just  been 
removed)  such  that  the  location  of  the  thermocouple  was  approximately  at  the  center  (i.e., 
along  the  axis)  of  the  sound  chamber.  This  thermocouple  temperature  was  then  monitored 
until  it  “plateaued”  and  provided  the  measurement  of  the  ambient  temperature  within  the 
sound  chamber. 

This  completed  a  single  experimental  trial  for  a  specified  resonant  frequency, 
pressure  ratio  and  temperature.  The  procedure  was  carried  out  a  total  of  nine  times  for  each 
different  pressure  ratio.  Once  the  six  outputs  for  each  run  were  recorded,  they  were  then 
transferred  over  to  a  spreadsheet  where  all  other  significant  parameters  were  computed 
automatically.  The  following  section  delineates  the  various  calculations  performed  in  the 
spreadsheet. 

E.  EXPERIMENTAL  CALCULATIONS 

The  first  calculation  desired  from  the  spreadsheet  entails  finding  the  actual  power 
being  supplied  to  the  cartridge  heater.  This  is  expressed  in  terms  of  the  current  and  voltage 
outputs  from  each  experimental  ran. 


P  =  IV 


(36) 
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Once  the  power  to  the  cartridge  heater  was  known,  the  surface  temperature  of  the  test 
cylinder  itself  could  be  calculated.  Knowing  the  resistance  of  the  thermal  circuit  as  given  in 
Appendix  B,  the  difference  between  the  interior  temperature  and  the  surface  temperature  is 
simply 


aT  =  PReq  (37) 

Utilizing  the  interior  temperature  output  as  provided  by  the  thermocouple  embedded  in 
the  cartridge  heater,  the  surface  temperature  of  the  cylinder  is  then  expressed  as 


T.  -  PR 

c  eq 


(38) 


Once  the  surface  temperature  is  known,  the  difference  between  it  and  the  ambient 
temperature,  as  measured  during  the  experiment,  is  calculated.  The  result  is  then  combined 
with  the  power  calculated  in  Eq.  35,  as  well  as  with  the  external  surface  area  of  the  test 
cylinder,  to  find  the  convective  heat  transfer  coefficient 


h  = 


P 

MTS  ~  Ta) 


The  Nusselt  number  is  then  derived  from  the  following  equation 


(39) 


Nu  =  —  (40) 

k 

where  d  is  the  test  cylinder  diameter  and  k  is  the  thermal  conductivity  of  air. 

The  next  step  is  to  calculate  the  various  criteria  as  previously  listed  in  the  theory 
section.  In  order  to  find  the  length  scale  ratio,  %,  the  speed  of  sound  within  the  chamber  must 
first  be  calculated  using  the  ambient  temperature 


c  =  JyR(Ta  +  273.15) 


(41) 
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The  value  of  x  can  now  be  found  from  Eqs.  1  -  4.  The  amplitude  parameter,  e,  is  derived  in 
the  spreadsheet  by  combining  the  pressure  ratio  from  Eq.  33  along  with  Eq.  8.  The  third 
parameter  calculated  is  the  frequency  parameter.  A2,  from  Eqs.  10  and  11. 

The  streaming  Reynolds  number,  Rs,  is  found  by  using  Eq.  16.  This,  however,  does 
not  represent  the  true  value  at  the  test  cylinder  position  due  to  its  not  being  precisely  at  the 
velocity  antinode.  A  corrected  value  for  Rs  can  be  deduced  by  finding  the  particle  velocity 
offset  between  the  velocity  antinode  location  and  the  test  cylinder  position.  This  offset  is 
derived  by  knowing  the  frequency  of  the  sinusoidal  signal  being  generated,  the  speed  of 
sound  within  the  chamber  from  Eq.  40,  and  the  distance  from  the  test  cylinder  to  the  end 
plate  termination  face,  L.  Hence 


t/0|sin(  —  -  )| 
01  c/f  2  1 


or 


Rs  =  ^Jsin(-—“)|2 

**  corrected 


from  which  it  was  found  that  the  corrected  value  for  Rs  had  less  than  a  0. 1%  error  due  to  the 
slightly  displaced  location  of  the  cylinder. 
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IV.  RESULTS  AND  DISCUSSION 


Nearly  600  experimental  trials  were  performed  throughout  the  course  of  this  study 
which  produced  results  for  183  distinct  data  points.  Data  were  obtained  for  five  different 
frequencies  at  various  pressure  ratios  ranging  from  0.9  %  to  3.2  %.  Three  separate  series  of 
trials  were  performed  at  each  acoustic  signal  setting,  (i.e.,  at  each  frequency  and  pressure 
ratio  setting)  corresponding  to  three  separate  settings  for  the  driving  temperature  difference 
between  the  test  cylinder  surface  and  ambient  conditions.  The  resultant  values  of  the 
streaming  Reynolds  number  ranged  from  40  to  1070  while  the  corresponding  Nusselt 
numbers  obtained  varied  from  8  to  38.  Figure  14  is  a  parameter  map  as  suggested  by 
Richardson  (1967)  which  shows  the  range  of  values  for  the  experimental  data  covered  plotted 
as  a  function  of  the  amplitude  parameter  versus  the  frequency  parameter,  and  delineates  the 
expected  different  regimes  of  flow.  The  data  obtained  cover  a  very  narrow  regime  of  this 
parameter  map  as  was  intentionally  planned  for  this  experiment.  Now  it  can  be  seen  that  the 
heat  transfer  results  obtained  through  experimentation  are  quite  evenly  distributed  between 
two  distinct  regions  on  the  map.  Region  A  represents  the  regime  in  which  the  flow  is 
expected  to  remain  laminar,  incompressible  and  attached,  and  outer  acoustic  streaming  is  the 
main  heat  transport  mechanism.  This  region  is  well  understood  in  theory  but  has  yet  to  be 
thoroughly  verified  through  experimentation.  It  is  anticipated  that  the  heat  transport 
mechanism  for  the  data  in  region  B  will  be  presented  as  a  combination  of  effects,  including 
that  of  vortex  shedding,  as  predicted  by  Honji  (1981),  Hall  (1984)  and  Sarpkaya  (1986). 

The  heat  transfer  results  are  presented  as  plots  of  Nu  vs.  Rs  (Figure  15).  Here  the 
difference  between  the  two  regions  of  varying  heat  transport  mechanisms  becomes  more 
apparent.  A  break  clearly  occurs  in  the  region  where  Rs  ~  240  (±  20)  and  it  can  be  concluded 
that  there  is  some  critical  point  in  this  range  where  there  is  a  transition  in  the  flow  at  which 
vortex  shedding  begins  to  become  a  dominant  factor  in  the  heat  transport  away  from  the  test 
cylinder.  In  order  to  better  examine  the  differences  between  these  two  regimes,  the  data  is 
divided  into  their  respective  groups  and  individually  analyzed. 
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Figure  14.  Parameter  Map  of  Expected  Heat  Transfer  Regimes  as  presented  by 
Richardson  (1963):  Convection  by  Inner  Acoustic  Streaming  (A),  by  Outer  Acoustic 
Streaming  (B). 
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Figure  15.  All  Data  Plotted  as  a  Function  of  Nusselt  Number  vs.  Streaming  Reynolds 
Number. 
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A.  LAMINAR,  ATTACHED  FLOW  REGIME 

Figure  16  is  a  plot  of  the  heat  transfer  results  in  terms  of  the  Nusselt  number  versus 
the  streaming  Reynolds  number  for  the  data  in  which  Rs  <  240,  and  for  which  criteria  A 
through  C  (Eqs.  4, 5, 15)  have  been  met.  Those  data  points  which  do  not  meet  these  criteria 
have  been  disregarded.  The  critical  parameters  for  the  remaining  data  are  as  follows: 

X<0.1 

e  <  0.3 

A2  >1800 

R 

-2  <  4.5 
A 

Theory  clearly  indicates  that  the  dependency  of  the  Nusselt  number  on  the  streaming 
Reynolds  number  in  this  regime  is  of  the  form  Nu  =  xPr^R,.05.  Since  the  Prandtl  number 
remains  constant  throughout  the  experiment,  the  solution  for  this  dependency  be  can  further 
simplified  as  Nu  =  CRS0  5,  where  the  term  “C”  encompasses  both  the  Prandtl  number  and  the 
qualitative  constant  of  the  previous  equation.  However,  this  solution  form  is  only  valid  for 
“large”  values  of  Rs.  Since  the  theory  does  not  provide  a  definite  limit  for  what  qualifies  as 
“large”,  this  criterion  had  to  be  determined  from  a  careful  examination  of  the  experimental 
data.  From  Figure  16,  it  can  be  observed  that  there  is  indeed  a  break  point  at  Rs  ~  130  where 
the  results  diverge  into  two  separate  solutions.  It  was  found  that  the  square-root  dependency 
on  Rs  does  not  significantly  change  past  this  value,  and  it  is  therefore  suggested  that  this  may 
be  in  the  range  of  the  lower  limit  of  “large”  values  of  Rs.  For  those  values  in  which  Rs>  130, 
a  curve  fit  of  the  heat  transfer  characteristics  results  in  a  solution  of  the  form 

Nu  =  0.94RS05  (45) 

and  it  can  therefore  be  determined  that  the  range  of  values  130  <  Rs  <  240  is  representative 
of  “large”  values  of  the  streaming  Reynolds  number.  The  values  below  Rs  =  130  are 
excluded  as  not  being  large  enough  due  to  a  variety  of  reasons  as  described  later. 
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Davidson  (1973),  in  an  extension  of  the  work  by  Richardson,  analytically  and 
numerically  tackled  the  problem  of  heat  transfer  from  a  cylinder  in  a  strong  acoustic  field  in 
great  detail.  He  obtained  a  correlation  of  the  dependency  of  the  Nusselt  number  on  both  the 
Prandtl  and  streaming  Reynolds  numbers.  The  correlation,  as  extracted  from  the  work  of 
Davidson  by  Gopinath  and  Mills  (1993)  for  this  regime,  is  of  the  form 

Nu  =  l388Pr0J%0-5  (46) 

By  taking  Pr  =  0.7  for  air  for  these  experiments,  the  equation  then  becomes 

Nu  =  1.07  R?-5  (47) 

The  experimental  fit  in  Eq.  45  under  predicts  by  about  13%,  but  supports  this  correlation 
well  within  the  limits  of  uncertainty. 

Figure  17  is  a  plot  of  the  region  in  which  Eq.  45  is  valid,  and  includes  the 
experimental  uncertainty  of  each  data  point  as  derived  in  Appendix  C.  It  can  be  observed 
from  this  plot  that  the  deviation  from  the  curve  fit  in  this  range  of  values  for  Rs  is  well  within 
experimental  uncertainty  limits. 

Although  the  lower  end  of  the  range  for  “large”  Rs  in  which  the  predicted  solution  is 
valid  has  been  determined  to  be  130  for  the  results  obtained  from  these  experiments,  it  is  by 
no  means  an  absolute  boundary.  Even  though  the  resultant  heat  transport  characteristics 
deviate  significantly  below  that  point  in  the  region  where  “intermediate”  values  of  R,  are 
present,  there  are  several  factors  which  could  account  for  part  of  the  discrepancy,  especially 
in  the  region  around  100  <  Rs  <  130.  These  include  uncertainty  due  to  equipment  limitations 
and  the  effects  of  natural  convection  and  conduction  of  heat  away  from  the  test  cylinder. 

The  effect  of  natural  convection  on  most  of  the  intermediate  values  of  Rj  is 
negligible,  though,  as  characterized  by  the  very  low  ratio  of  the  Grashof  number  to  the  square 
of  the  streaming  Reynolds  number  (except  at  veiy  low  values  of  Rs).  The  effect  due  to 
conduction  is  much  harder  to  quantify  in  so  simple  a  form,  though.  Conduction  can  and  does 
occur  during  the  experiment  at  two  separate  places  where  the  test  cylinder  is  in  contact  with 
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other  apparatus.  The  top  of  the  test  cylinder  is  in  contact  with  the  “plug”  used  to  hold  it  in 
place,  while  the  bottom  of  the  test  cylinder  is  allowed  to  rest  in  a  small  depression  on  the 
wall  of  the  plexiglass  sound  chamber.  Realizing,  though,  that  the  contact  area  at  both  points 
is  very  small,  and  that  the  thermal  resistance  of  the  plug  and  the  plexiglass  wall  are  both 
relatively  high  (due  to  their  low  thermal  conductivity),  the  heat  transport  away  from  the  test 
cylinder  can  be  assumed  to  be  negligible  as  well. 

The  uncertainty  in  the  results  due  to  the  equipment  limitations,  though,  does  have  a 
very  significant  impact  on  the  results.  The  calculated  value  of  the  Nusselt  number  has  an 
uncertainty  of  up  to  20%  (  see  Appendix  C)  depending  upon  the  power  being  dissipated  by 
the  test  cylinder.  Since  the  experiment  revolves  around  finding  the  heat  transfer 
characteristics  at  specific  values  of  temperature  difference  between  the  test  cylinder  surface 
and  the  ambient  conditions,  correspondingly  low  power  dissipation  from  the  cylinder  occurs 
as  the  streaming  Reynolds  number  decreases.  Therefore,  the  region  of  intermediate  values 
of  Rs  have  relatively  low  electrical  heat  dissipation,  and  hence  low  current  values  associated 
with  them,  dropping  to  as  low  as  0.06  amps  in  some  cases.  Since  the  equipment  uncertainty 
for  the  current  reading  is  of  the  order  of  the  last  digit  present,  this  particular  component  of 
the  Nusselt  number  has  an  uncertainty  of  nearly  18%  by  itself  and  greatly  influences  the 
overall  uncertainty.  Therefore,  it  may  be  more  accurate  to  define  the  lower  limit  of  large  Rs 
values  as  somewhere  between  100  and  130.  However,  the  error  due  to  equipment 
measurements  is  not  enough  to  compensate  for  the  disparity  between  theory  and  experiment 
at  values  much  less  than  100. 

B.  SEPARATED  FLOW  REGIME 

The  second  regime  which  this  experiment  encompasses  is  that  in  which  vortex 
shedding  and  other  forms  of  unsteady  flow  begin  to  affect  the  heat  transport  characteristics. 
Figure  18  is  a  plot  of  the  resultant  data  in  this  regime  as  obtained  during  experimentation  in 
terms  of  the  Nusselt  number  versus  the  streaming  Reynolds  number.  A  curve  fit  of  the  data 
results  in  a  solution  of  the  form 


Nu  =  0.31RS069 


(48) 
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Figure  18.  Unstable  Regime. 
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There  is  no  theory  for  this  flow  regime  with  which  to  compare  the  results,  but  it  is  reasonable 
to  assume  that  this  curve  fit  is  representative  of  the  correct  solution.  As  expected,  the 
unstable  flow  resulting  at  higher  values  of  Rs  would  increase  the  heat  transfer  rate  from  the 
cylinder,  and  hence  the  stronger  dependency  that  the  Nusselt  number  has  on  the  streaming 
Reynolds  number,  as  opposed  to  the  usual  square  root  dependence. 
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V.  CONCLUSIONS 


Experiments  were  conducted  to  observe  the  convective  heat  transfer  rates  to  an 
isolated  cylinder  in  an  acoustic  standing  wave.  A  comparison  of  the  various  length  scales 
and  other  parameters  was  conducted,  and  the  experimental  method  stated.  During  the 
experiment,  the  properties  of  the  acoustic  field  were  varied  to  provide  a  large  base  of  data 
which  was  then  analyzed  and  discussed.  Several  regimes  of  interest  were  investigated  and 
the  results  presented.  Figure  19  is  a  plot  of  all  data  obtained  with  curve  fits  for  the  regimes 
of  interest. 

Essentially,  heat  transfer  from  a  cylinder  in  a  zero-mean  oscillatory  flow  as 
represented  by  an  acoustic  standing  wave  can  be  divided  into  at  four  separate  regimes  in 
which  different  heat  transport  mechanisms  dominate.  For  very  low  values  of  the  streaming 
Reynolds  number,  Rs,  convective  effects  due  to  the  acoustic  field  are  negligible  and  natural 
convection  is  then  the  dominant  mode  of  heat  transport.  For  intermediate  values  of  Rs,  there 
is  a  stronger  dependence  on  Rs  but  not  yet  on  the  order  of  Rs  0  5  since  Rs  is  still  not  large 
enough  for  flow  to  be  of  the  boundary  layer  type.  Buoyancy  effects  are  comparable  in  the 
lower  end  of  this  regime,  becoming  small  for  larger  values  of  Rs.  For  much  larger  values  of 
Rs,  past  100  or  so,  an  acoustic  streaming  flow  presents  itself  in  the  boundary  layer,  resulting 
in  a  square-root  dependency  on  Rs.  Experimentally,  the  results  obtained  in  this  regime 
closely  match  the  expected  theory,  and  the  heat  transfer  characteristics  may  be  estimated  by 
Eq.  45.  Finally,  past  a  critical  value  of  Rs  -240  (which  confirms  well  with  theory),  an 
unstable  flow  with  vortex  shedding  begins  to  take  place  at  the  surface  of  the  cylinder, 
increasing  the  dependency  on  Rs  which  the  overall  heat  transfer  solution  has.  The  heat 
transfer  characteristics  in  this  regime  may  be  estimated  from  Eq.  48.  It  is  these  last  two  of 
the  above  regimes  that  formed  the  focus  of  this  study. 
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VI.  RECOMMENDATIONS 


Several  directions  for  further  study  may  be  suggested  by  the  data  obtained  from  this 
work.  One  factor  which  may  need  further  investigation  is  a  study  of  the  effects  of  varying 
the  aspect  ratio  of  the  chamber  diameter  to  the  cylinder  length  to  determine  whether  it  was 
large  enough  for  these  experiments.  A  ratio  of  15  was  used  for  this  experiment,  and  was 
assumed  to  be  large  enough  to  discount  flow  effects  caused  by  the  walls  of  the  chamber,  but 
this  number  was  chosen  somewhat  arbitrarily  and  only  a  detailed  experimental  study  can 
determine  the  actual  effect. 

Another  possible  source  of  error  which  cannot  be  accurately  accounted  for  concerns 
itself  with  the  experimental  method  used.  Specifically,  the  measurement  of  the  ambient 
temperature  within  the  chamber  at  the  time  of  each  experimental  trial.  The  method  used,  that 
of  removing  the  test  cylinder  from  the  sound  chamber  after  its  temperature  has  been  recorded 
and  the  acoustic  signal  has  ceased,  then  placing  a  thermocouple  through  the  hole  it  has  just 
evacuated  in  the  chamber  wall,  can  be  improved  upon.  During  the  time  required  while 
waiting  for  the  thermocouple  temperature  to  peak,  the  heated  air  in  the  sound  chamber  is 
rising  to  the  top  of  the  chamber  and  the  accuracy  of  correlating  the  thermocouple  temperature 
reading  to  that  of  the  ambient  temperature  at  the  time  of  the  experimental  trial  is  left  in 
doubt.  A  better  method  would  be  to  permanently  affix  the  thermocouple  in  the  chamber  so 
that  simultaneous  measurement  of  both  the  cylinder  temperature  and  the  ambient  temperature 
can  be  taken. 

Additional  research  into  three  different  areas  of  the  problem  come  to  mind.  First,  a 
test  of  the  effects  of  placing  the  test  cylinder  horizontally  in  the  sound  chamber  is  suggested, 
although  there  should  be  little  if  any  dependency  on  this  orientation  since  natural  convection 
effects  are  small  for  the  strong  acoustic  fields  being  used.  Completely  new  geometries  may 
also  be  tested  which  would  mimic  actual  heat  exchanger  component  shapes  expected  in  a 
thermoacoustic  engine.  Finally,  additional  research  using  different  gases  in  the  sound 
chamber  would  provide  data  on  the  dependency  of  the  Nusselt  number  on  the  Prandtl 
number. 
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APPENDIX  A.  CALIBRATIONS  AND  CALCULATIONS 


Several  pieces  of  the  experimental  apparatus  required  some  form  of  calibration,  or 
calculation  of  a  specific  parameter,  prior  to  initiating  the  experiments.  The  most  important 
equipment  items  of  concern  were  the  “unattached”  J-type  thermocouples  which  were  to  be 
used  for  various  applications  throughout  the  experiment,  as  well  as  the  thermocouple  which 
was  embedded  in  the  cartridge  heater.  In  addition,  an  equivalent  thermal  resistance  for  the 
cartridge  heater/silicon  oil/copper  sheath  circuit  needed  to  be  calculated  along  with 
assurances  that  the  linear  temperature  distribution  along  the  test  cylinder  was  within 
reasonable  limits.  An  additional  study  was  performed  to  analyze  the  heat  transfer  effects  on 
the  test  cylinder  due  only  to  natural  convection. 

Three  J-type  thermocouples  were  used  throughout  the  experiment  and  were  the  first 
items  to  be  calibrated.  Since  accurate  temperature  information  was  crucial  to  the  reliability 
of  the  data  obtained  through  experimentation,  the  thermocouples  were  tested  to  see  if  any  of 
them  showed  a  tendency  to  read  either  higher  or  lower  than  the  actual  temperature  (a 
somewhat  common  occurrence).  The  embedded  thermocouple  in  the  cartridge  heater  was 
also  tested  for  the  same  reason.  All  thermocouple  leads  were  attached  to  the  thermocouple 
reader  to  be  used  throughout  the  experiment  to  ensure  that  the  entire  circuit  was  tested 
concurrently. 

The  unattached  thermocouples  and  the  cartridge  heater  were  all  placed  in  an  ethyl 
glycol  solution  belonging  to  a  Rosemont  Model  913A  calibration  bath.  A  Rosemont  Model 
920A  commutating  bridge,  which  utilizes  a  precision  temperature  probe  as  its  input,  provides 
the  reference  temperature.  The  temperature  of  the  bath  was  then  set  at  varying  points 
between  22°  and  48°C,  the  expected  experimental  temperatures  being  well  within  that  range. 
Table  2  shows  the  results  of  the  calibration.  All  thermocouples  were  found  to  read  within 
0.1  °C  of  the  reference  temperature  for  all  cases.  This  was  well  within  the  possible 
uncertainty  of  0.5°C  listed  for  J-type  thermocouples. 

Next,  Figure  20  shows  the  thermocouple  arrangement  which  was  used  to 
experimentally  derive  the  equivalent  thermal  resistance  and  the  linear  temperature 
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Reference 

Cartridge 

Thermocouple 

Thermocouple 

Thermocouple 

Maximum 

Temperature 

Heater 

“A” 

“B” 

“C” 

Deviation 

(C) 

(C) 

(C) 

(C) 

(C) 

(C) 

22.83 

22.8 

22.8 

22.8 

22.8 

<0.1 

23.96 

23.9 

23.9 

23.9 

23.9 

<0.1 

26.83 

26.8 

26.8 

26.8 

26.8 

<0.1 

29.63 

29.6 

29.6 

29.6 

29.6 

<0.1 

32.49 

32.5 

32.5 

32.5 

32.5 

<0.1 

35.37 

35.4 

35.4 

35.4 

35.4 

<0.1 

38.04 

38.1 

38.0 

38.0 

38.0 

<0.1 

40.74 

40.8 

40.8 

40.8 

40.8 

<0.1 

43.02 

43.0 

43.0 

43.0 

43.0 

<0.1 

45.68 

45.7 

45.7 

45.7 

45.7 

<0.1 

48.19 

48.2 

48.2 

48.2 

48.2 

<0.1 

Table  2:  Thermocouple  Calibration  Data  (all  values  in  °C) 


Figure  20.  Thermocouple  Arrangement  for  Calibration  Tests. 
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distribution  data.  The  three  unattached  thermocouples  which  were  previously  calibrated 
were  securely  placed  on  the  outer  surface  of  the  test  cylinder  using  clamps,  the  test  cylinder 
being  already  prepared  as  it  normally  would  for  an  experimental  run.  The  test  cylinder  was 
then  suspended  horizontally  and  power  was  supplied  to  the  cartridge  heater.  Once  all  four 
of  the  thermocouples  reached  a  steady-state  temperature,  the  temperatures  were  recorded 
along  with  the  voltage  and  current  being  supplied  to  the  heater.  The  equivalent  thermal 
resistance  was  then  derived  using  the  equation 


R 


eq 


T  -  T 

_ C _ S 

IV 


(A.l) 


Temperatures  which  were  to  be  comparable  with  those  used  during  the  experiments  were 
obtained  in  a  still  air  environment.  Some  additional  tests  were  performed  using  a  fan 
blowing  air  across  the  test  cylinder  to  create  a  higher  heat  transfer  rate  away  from  the 
cylinder  so  that  higher  power  levels  to  the  heater  could  be  reached  while  staying  within  the 
limits  of  the  expected  temperature  range.  An  average  temperature  was  then  obtained  using 
all  three  thermocouples,  and  this  value  was  in  turn  utilized  to  derive  an  equivalent  thermal 
resistance  for  each  test  run.  The  final  value  of  the  thermal  resistance  resulted  in  averaging 
the  values  from  each  test  run,  obtaining  1 .022  K/W  with  a  maximum  deviation  of  0.089 
K/W.  The  results  are  annotated  in  Table  3. 

The  same  test  results  were  used  to  determine  the  change  in  temperature  along  the 
length  of  the  test  cylinder.  The  difference  in  temperature  along  the  length  at  the  lower  power 
settings  demonstrated  a  fairly  even  distribution  of  heat,  with  a  variation  of  less  than  0.6°C. 
This  variation  increased  as  the  power  to  the  heater  was  also  increased,  as  could  be  expected. 
For  the  temperature  range  in  which  the  experiments  were  run  (again,  26°  to  40°  C),  a 
maximum  variation  of  1.6°C  occurred,  although  this  was  at  a  very  high  level  of  power  being 
supplied  to  the  heater  and  was  not  representative  of  the  remaining  data. 
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Current 

(Amps) 

Voltage 

(Volts) 

Power 

(Watts) 

Interior 

Temp 

(C) 

Thermo¬ 
couple  “A” 

(C) 

Thermo¬ 
couple  “B” 

(C) 

Thermo¬ 
couple  “C" 

(C) 

Average 

Surf  Temp 

(C) 

Equiv 

Resistance 

(C) 

0.08 

8.30 

0.664 

38.4 

37.8 

38.1 

37.3 

37.73 

1.009 

0.07 

7.40 

0.518 

36.0 

35.5 

35.5 

35.4 

35.47 

1.023 

0.08 

8.40 

0.672 

37.3 

36.9 

36.8 

36.3 

36.67 

0.938 

0.09 

9.50 

0.855 

40.8 

40.2 

40.1 

39.5 

39.93 

1.018 

0.09 

10.30 

0.927 

42.9 

42.3 

42.0 

41.5 

41.93 

1.046 

0.10 

10.70 

1.070 

44.8 

44.2 

43.7 

43.2 

43.70 

1.028 

0.08 

9.00 

0.720 

27.4 

26.9 

26.6 

26.3 

26.60 

1.111 

0.09 

10.00 

0.900 

28.6 

27.9 

27.8 

27.2 

27.63 

1.078 

0.10 

11.00 

1.100 

30.1 

29.3 

29.2 

28.6 

29.03 

0.973 

0.11 

12.00 

1.320 

31.5 

30.5 

30.5 

29.6 

30.20 

0.985 

0.14 

15.50 

2.170 

37.6 

35.7 

36.0 

34.4 

35.37 

1.028 

Table  3:  Thermal  Resistance  and  Linear  Temperature  Distribution  Trials 
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APPENDIX  B.  UNCERTAINTY  ERROR  ANALYSIS 


In  order  to  obtain  a  measure  of  the  reliability  of  the  data  obtained  during 
experimentation,  an  uncertainty  error  analysis  was  performed  for  both  the  Nusselt  number 
and  the  streaming  Reynolds  number.  The  analysis  consisted  of  finding  the  maximum 
possible  deviation  for  all  components  appearing  in  the  equations  which  define  both 
parameters.  Then,  in  a  standard  fashion,  a  root  mean  square  analysis  was  performed  to 
derive  a  reasonable  overall  possible  error  for  each  experimental  run.  The  error  analysis 
formulae  were  themselves  incorporated  into  the  “results  worksheet”  provided  in  Appendix 
C  so  that  each  data  point  has  an  associated  possible  error  derived  from  the  input  provided. 

The  Nusselt  number  can  be  derived  as  shown  from  Eqs.  37  through  40 


Nu 


_ IV _ 

nkl(T  -  IVR  -  T) 

v  c  eq  a ' 


(B.l) 


where  l  is  the  length  of  the  test  cylinder. 

Using  the  manufacturers’  recommended  equipment  error  ranges,  the  maximum 
uncertainty  in  each  measured  component  in  Eq.  B.l  is  as  follows 


V  =  0.05%  reading  +  0.02%  full  scale  +  1  digit 

(B.2) 

I  =  0.5%  RDG  +  1  digit 

(B.3) 

Tc  =  ±0.5"  C 

(B.4) 

Ta  =  ±0.5"  C 

(B.5) 

The  maximum  uncertainty  in  the  equivalent  resistance  is 

calibration  data  in  Appendix  A  as 

obtained  from  the 
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R  =  ±0.089  K/W 

eq 


(B-6) 


The  analysis  for  the  overall  uncertainty  itself  is  structured  as  follows.  For  the  Nusselt 
number,  the  root  mean  square  error  is  given  by 


a  m  =  [E 

1  OX; 


(B-7) 


where  X;  represents  each  individual  component  of  Eq.  B.l.  The - term,  as  the  partial 

dX. 

derivative  indicates,  physically  represents  the  sensitivity  of  the  Nusselt  number  to  the 
variable  Xj,  provided  all  other  variables  are  unchanged.  The  aX.  represents  the  uncertainty 
in  the  corresponding  variable  as  given  in  Eqs.  B.2  to  B.5.  For  instance,  the  contribution  to 
the  uncertainty  due  to  the  voltage  measurement  is 


(A  Nu)y 


=  Nu - — - Ta)(A^ — 

av  V(Tc  -  VIReq  -  Ta) 


(B-8) 


A  value  with  more  significance,  though,  is  the  individual  fractional  uncertainty  which 
takes  into  account  the  calculated  value  of  the  Nusselt  number  and  can  be  expressed  as  a 
percentage  possible  error.  This  is  represented  by  dividing  the  individual  uncertainty  by  the 
Nusselt  number  in  the  following  manner 


.  (aM0„  (Tc  ~  Ta)( aV) 

individual  fractional  uncertainty  =  -  =  - 

Nu  V(Tc  -  VIReq  -  Ta) 


(B.9) 


This  new  term  leads  directly  to  the  desired  method  of  expressing  the  possible  error 
in  the  calculated  value  of  the  Nusselt  number  as  an  overall  fractional  uncertainty  using  a  root 
mean  square  analysis. 


A  (NU)  _  r£(^)M 
Nu  i  Nu 


(B.10) 


In  a  similar  fashion,  the  individual  fractional  uncertainty  of  the  remaining  terms  are 
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as  follows 


(aNu)j 

c Te 

-  Ta){Al) 

Nu 

I(TC  - 

W*"  ~ 

T<) 

(a  Nu)t 

W 

Nu 

(L  - 

-  WR*  - 

Ta) 

(aNu)t 

(a  T) 

v  a' 

Nu 

or  - 

-  WR.,  ~ 

Ta) 

(  aNu)^ 

Nu 

Vc  - 

-  WR*  - 

Ta) 

(B  .11) 

(B.12) 

(B.13) 

(B.14) 


By  examining  the  data  results,  the  largest  contributor  to  the  overall  error  in  the  Nusselt 
number  is  due  to  the  current.  This  is  caused  by  the  small  currents  being  utilized  during  the 
experiment,  which  were  as  low  as  0.05  amps.  Since  the  uncertainty  in  the  current  is  of  the 
order  of  0.01  amps,  there  can  be  an  error  in  the  calculated  Nusselt  number  of  approximately 
20%  due  to  the  current  term  alone.  One  way  to  lessen  the  effect  that  the  current  term  has  on 
the  overall  error  is  by  decreasing  the  voltage  output  of  the  power  supply  and  hence  increasing 
the  current  needed  to  maintain  the  same  power  being  generated.  This  capability,  though,  is 
not  a  feature  of  the  equipment  being  used.  It  must  also  be  noted  that  although 
a Tca  =  ±0.5°  Cin  Eqs.  B.4  and  B.5,  the  calibration  of  the  thermocouples  described  in 
Appendix  A  indicated  an  error  of  less  than  0. 1°C  for  the  temperature  range  of  the  experiment 
If  this  is  taken  into  account,  the  error  due  to  the  ambient  temperature  (Ta)  and  the  center 
temperature  (Tc)  terms  in  Eq.  B.l  would  diminish  by  80%. 

The  error  analysis  for  the  streaming  Reynolds  number  is  similar  to  that  just  performed 
for  the  Nusselt  number.  Utilizing  Eq.  17  and  substituting  Eqs.  33  and  41  into  it,  the  value 
of  Rs  can  be  shown  to  be 


R(Ta  *  273.15)(y0V2)2 
2Tt/vyPm2(0.523)2 


(B.15) 


Recall  that  Vjfe  is  the  multimeter  output  voltage  as  derived  from  the  pressure  transducer 


after  being  passed  through  the  100  gain  preamplifier.  Thus,  Rs  becomes  a  function  of  only 
the  following  measured  variables:  the  frequency,  the  ambient  temperature  and  the  pressure 
transducer  output  as  read  on  the  multimeter,  while  all  additional  parameters  remain  constant. 
The  maximum  uncertainty  for  the  thermocouple  is  the  same  as  previously  listed  (±0.5°C) 
while  the  other  two  variables  have  uncertainties  provided  in  the  manufacturers’ 
specifications.  Calibration  data  for  the  pressure  transducer  used  during  the  experiment 
indicates  an  uncertainty  equivalent  to  0.15%  of  the  percentage  of  Full  Scale  Output  (FSO) 
where  the  FSO  is  254  mV.  This  correlates  to  a  maximum  deviation  of  2.3  mV  after  passing 
through  the  100  gain  preamplifier.  In  addition  to  this,  the  multimeter  which  is  used  to  read 
this  voltage  has  an  error  of  0.06%  reading  +  0.03%  Range.  These  combine  to  give  a  total 
microphone  voltage  output  error  of 

Mic  =  2.3  mV  +  0.06%  reading  +  0.03%  range  (B.16) 

The  uncertainty  in  the  frequency  signal  from  the  function  generator  error  is  given  as  20  ppm, 
i.e., 

/  =  20x10 -6  reading  (B.17) 


Again  a  root  mean  square  analysis  was  performed  in  a  manner  similar  to  Eq.  B.7  to  derive 
the  overall  fractional  uncertainty  for  the  streaming  Reynolds  number. 


_  r£ 

R  i  R 


(B.18) 


The  individual  fractional  uncertainties  are 

_  a?; 

Rs  ~  (Ta  +  273.15) 
(a  RJmc  =  2(a  Mic) 

Rs  yV2 


(B.19) 

(B.20) 
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(B.21) 


Af 

^  / 

The  largest  contribution  to  the  error  in  the  streaming  Reynolds  number  is  the  error 
due  to  the  microphone  output  voltage,  specifically,  the  possible  error  in  the  pressure 
transducer  itself.  However,  the  data  presented  in  Appendix  C  shows  that  this  error  is  very 
small  and  in  the  range  of  <2%  of  the  total  value. 
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APPENDIX  C.  EXPERIMENTAL  DATA 


The  following  pages  contain  the  data  obtained  through  experimentation  in 
spreadsheet  fashion.  All  of  the  relevant  parameters  are  listed,  although  some  constants  have 
been  left  out  due  to  size  constraints. 
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